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INVERSION OF SOME LAPLACE TRANSFORMS 
ENCOUNTERED IN HEAT TRANSFER PROBLEMS. 

by 

JOHN E. BROCK 

Professor of Mechanical Engineering, NPS 
ABDOLLAH ZARGARY 

Lieutenant, Imperial Iranian Navy 
and 

CRAIG COMSTOCK 

Associate Professor of Mathematics, NPS 
Introduction 

The Laplace inversions treated herein arise in connection with some 
conjugated heat transfer problems, details of which have been reported in 
a thesis by one of the authors [4] and will also be reported elsewhere, 
as well. In connection with another physical problem, R. A. Kantola [i] 
arrived at similar results. However, the details here are sufficiently 
involved and sufficiently different than those dealt with by Kantola to 
suggest preserving our analysis in a monograph such as this. 

In the body of this monograph, a theorem is stated and proved which 
affords the Laplace inversion of a rather general function. Appendix A 
indicates a suggested numerical procedure for evaluation of the definite 
integrals appearing in the result. Appendix B treats certain details of 
a special but useful case. Appendix C merely lists some additional Laplace 
inversions which have come to our attention in connection with related work 
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but which do not result from the analysis given in the body hereof. 



Statement of Theorem 

The principal content of this monograph is contained in the 
THEOREM : Let 

f(s) = 7 exp[-/ tn(s) ] 



( 1 ) 



where m(s) is a meromorphic function whose zeros s = -aj^ and whose poles 



s = -b|^ satisfy 



0 < a-| < b-| < a2 < b^ < . . . < a|^ < b|^ i 



( 2 ) 



^ 1 

Then, for t > o, the inverse Laplace transform F(t) =«0 {f(s)} 
is given by N bj^ 

F(t) = exp [-y(o)] ■ ^ ^ / X [u(x)] dx (3) 



where 



k=l a, 



y(x) = /|m(-x) I 



(4) 



The sum, having upper limit N, is to be interpreted as follows. 

If there is a finite number N of pairs aj^ < bj^, no special interpretation 
is necessary. If there is a finite number (N-1) of pairs aj^ < b|^ and an 
additional zero s = -a|^, aj^| > , an Nth pole s = -b^^^ = - °° is to be 

added so as to make N pairs. If there is an infinite number of pairs 
aj^ < b, take N = «>. 



Proof of Theorem 

In the proof we consider a sequence {R^} of positive numbers, 

Rn ^ ^k’ *^k ^ R^ as n and we consider a 

corresponding sequence of closed contours in the complex s-plane each 
of which encloses the pole of f(s) at s = 0 but excludes all of the branch 
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points at s = -a 1^ 

contours C^-j , C^2» 

s = -i R to s = i 
n 



and s = -b|^. The contour consists of four sub- 
C and C . as follows. C , is a vertical line from 



R^, indented as shown by a small semicircle passi 



ng 




to the right of s = 0 . C^2 ^ |- = e ^ ir and is 
the similar arc s = -tt = 0 = - Finally, contour 
which passes to the right from s = -R^, along the negative real s-axis 
to s = -a.| with small semicircular indentations into the upper half 
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s-plane about each zero (-a|^) and about each pole (-bj^) after which it 

then passes to the left from s = -a, to s = -R with small semicircular 

I n 

indentations into the lower half s-plane, as indicated in Figure 1. 
Clearly, by Cauchy's theorem 

1 r 






e^^ f(s) ds = exp [-/ m(0) ] = exp [-m(o)] (5) 



2-ni 



J 



for n = 1 , 2, 3, ... Also, by the complex inversion theorem for Laplace 
transforms ^ 



1 



nl 27ri 



e^^ f(s) ds ^ F(t) 



( 6 ) 



j 



-nl 



as n ^ 



To proceed further, we wish to establish an inequality for arg[m(s)]. 
To do this, we consider complex vectors emanating from the points s = -aj^ 
and s = -b|^, and we write 



where 



and 



1 + s/a^ = U|^ exp(i4>,^); u,^ > 0‘. -ir < (|),^ = t 


(7) 


1 + s/b|^ = V|^ exp(it^i^); v,^ | q; -^ = ^ 


(8) 


(s) = m(o) (1 + s/a,^)T (1 + s/b,^) 


(9) 


2 

= (s) exp (i6n) 


(10) 


ul (s) = y(0) u,^ V v^ 


(11) 




(12) 



The sums and products written with a superscript (n) indicate that 
only those terms are to be considered for which a^ < or b^ < R^, as may 
be appropriate. 
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Now note for any s in the (closed) upper half plane we have 



'^k* l^“l> 2, ... 

A geometrical argument is given in Figure 2. Now either there are N 

pairs aj^, b|^ each less than R^, so that 

N 

k=l 

or there is an extra, unpaired zero, so that 

N-1 



( 13 ) 



(14) 



^n = 



(<i>k ■ = 0 



(15) 




In the first case, we can write 

N-1 

^n ^ ‘*’1 ■ k “ '*’k+l ) ■ 'I'n " "^l ■ < <l>i = TT 

k=l 

whereas in the second case we can write 

N 

®n " ‘*’1 ■ Z ^'•'k ■ “^k+l ) = 

Thus, in either case. 



k=l 



0 <= 6 , <= - 



(16) 



(17) 



(18) 



Similarly, by considering an s in the lower half plane, we find 



^ 0 



(19) 
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so that in all cases 



-TT ^ 1 IT 

Thus, if we write 

J - -L 

'^n2 " 2Tri 

we note that 

I f(s) I = ^ I exp (- / Pn2gi8 ) | 
n 

= ^ exp [-u^ cos(3^/2)] | ^ 




( 20 ) 



( 21 ) 



( 22 ) 



uniformly on virtue of the inequality (20). Thus, by Jordan's 

lemma 0^2 rn^. Similarly, the integral on vanishes as n ^ 



n4 



Next we consider the integrals around semicircles surrounding 
poles and zeros on contour The result is obvious near a zero: 
s = -a. ; e^^ e'^k^, ] f(s) | - since m(s) == 0» and the length 



of the arc goes to zero as its radius does. Thus, the integral on 
such a semicircle about a zero, s = -aj^, vanishes with the radius. 
However, for the integral on a semicircle about a pole, s = -b|^, it 
becomes necessary to use a considerably more subtle argument. For 
this purpose we must invoke (20) again. We have 



I exp[-m(s)] I = exp[-M^ cos(3^/2)] < 1 (23) 

even though -> «> as the radius of the indentation goes to zero. 

Thus, the integral on such a semicircle about a pole, s = -bj^, 
vanishes with the radius. 
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We remark that although we have reason to believe, on physical 
grounds, the correctness of an evaluation by Kantola [l], we feel 
that his argument intended to show vanishing of such an integral on a 
semicircle surrounding a branch point is deficient. Our argument, 
based on (20), is not applicable to his case. 

If now we write (5) in the form 



*^01 *^n ^^02 '^nS ^^n4 



(24) 



and let n and also let the radii of all semicircular indenta- 
tions on 9° to zero, we find 



F(t) = exp[-u(0)] - (25) 

where 

J - - 1 - 

*^3 ~ 2Tri 

and now is a contour extending from s = -“ along the upper edge of 
a branch cut, to the right, passing around s = -a-j, and then extending 
along the lower edge of the branch cut to s = The branch cut is 
the negative real s-axis from to -a^ . There are no longer any semi- 
circular indentations in C^. 

To reduce (26) to a more convenient form, we first introduce 
the function v(x), for real x, which denotes the excess of the number 
of zeros (s = -aj^) over the number of poles (s = -b|^) to the right of 
the point s = X. It is clear that v = 0 in the interval < s < -bj^ 

while V = 1 in the interval -bj^ < s < -aj^. Writing g now, rather than 
since we have taken n -> «> and all zeros and poles are included in 
our accounting, we note that on the upper edge of the cut B = vtt/2 



^ exp[st-/m(s) 



] ds 



(26) 
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while on the lower edge 6 =-vtt/2. We also now write y rather than y . 

n 

Thus, 



1 1 

3 “ 2ttT J - y exp(iv7T /2)]ds 



Jo = 



Ziri 



J. 



j exp[st - y exp(-ivy/2)]ds 



■ai 



-an 



27ri 



1 „st 

s ® 



exp[-y exp(ivir/2)] -exp[-y exp(-ivTr/2)] 



( 27 ) 



ds 



If V = 0, the quantity in braces vanishes, while if v = 1, it becomes 



exp(-iy) - exp(iy) = -2i sin 



(28) 



Thus 



■5 IT 



J e^^ si n y ds 



k=l 



If we now change the variable of integration from s to x = -s, 
we obtain 



J, = i 



N 



k=l 



^ e“^^ sin[y(x)]dx 



and, thus, finally 



F(t) = exp[-y(0)] - - 



' r 

-I / 

kr \ 



^ e'^* sin[y(x)]dx 



(29) 



(30) 



The use of the upper limit of summation is consistent with the remarks 
following the statement of the theorem, the proof of which is now complete. 
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APPENDIX A 



Numerical Evaluation of the Integrals 
The integrals appearing in the inversion theorem (3) are of the 



form 



where 



I = h(x)sin[g(x)] dx 
a 



(A-1) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 



h(x) is continuous in [a,b] 
g(x) is monotonic continuous in [a,b) 
g(a) = 0 

g(x) ^ as X -> b" 

It is contemplated that b may be finite or infinite. 

The integrand is illustrated in 
the figure at the right. The integral 
I, is the algebraic sum of the shaded 
areas, with signs taken as indicated. 

For values of x sufficiently close to 
b, these contributions will alternate 
in sign so that I may be regarded as 
the sum of an alternating series of 
terms having decreasing magnitude. 

Thus, the error of truncation is less than the magnitude of the first term 
to be neglected. 

In the numerical treatment one must evaluate, in succession, the 




numbers such that 



g (C^) = nTT 



The values evidently satisfy 



a = Cq < C-| < C 2 < . . . <b 



(A-6) 



(A-7) 
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This determination is uniquely possible because the monotonicity 
of g(x) assures the existence of a unique inverse function. Many tech- 
niques are available for this evaluation. Some remarks suggesting a 
procedure suited to the problems which form the subject of this mono- 
graph will appear near the end of the present appendix. 

Thus, we have 



where 




qj^ = I h(x) sin[g{x)] dx 




(A-8) 



(A-9) 



Numerical quadrature is indicated. A formula of Gaussian type 
might be appropriate except that the ordinary Gaussian points are all 
interior to the interval of integration and thus cannot reflect our 
knowledge that sin[g(C^)] = 0 for n = 0,1,2,... A slightly more ef- 
ficient procedure, accordingly, employs Gauss-Radau integration in 
which the end points of the interval are specifically included among 
the integration points. The integrand is zero at these points so that 
the corresponding summands may be omitted. For details of Gauss-Radau 
formulas, see Z. Kopal [2]. The appropriate formula is 



N-1 



F(x)dx = H,F(-1) + ^ F(aj) + F(l) (A-10) 



where a . = and Hj 



J=Z 

sre tabulated in Appendix 4.5 of 



Kopal [2] for N = 2,3,...,11. (Note that these a.'s should not be con- 

J 

fused with aj^'s which locate the zeros of the function m(s) employed in 
the body of this monograph). 
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For example, for n = 5, we have 



GAUSS-RADAU CONSTANTS FOR N = 5 
ag = 1 = -a^ ^5 = = ’’/''O 

a^ = / 3/7 = -a 2 = H 2 = 49/90 

83 = 0 H 3 = 64/90 

We have calculated the coefficients a. and H. for the case N = 12 

J J 

and list the values here so as to supplement the information given in [2]. 



GAUSS-RADAU CONSTANTS 


FOR N = 12 








a ^2 = 1 .00000 


00000 


00000 


0 


Hi 2 


= 


0.01515 


15151 


51515 


2 


a^^ = 0.94489 


92722 


22882 


0 


^11 


= 


0.09168 


45174 


13204 


7 


a^Q = 0.81927 


93216 


44007 


4 


*^10 


= 


0.15797 


47055 


64367 


1 


a^ = 0.63287 


61530 


31861 


8 


H9 


= 


0.21250 84177 


61024 


6 


ag = 0.39953 


09409 


65348 


9 


^^8 


= 


0.251 27 


56031 


99201 


1 


a-j = 0.13655 


29328 


54927 


6 


«7 


= 


0.27140 


52409 


10696 


2 



These formulas lead to exact evaluation for a polynomial F(x) of 
degree 2N-3. In our case, with appropriate mapping of [C^, onto 

[-1 , 1], we have F(-l) = F(l) = 0 so that using N = 9, say, requires 
seven nontrivial evaluations, leading to exact integration of a 15-degree 
polynomial, whereas with ordinary Gaussian integration, we would have had 
to use N = 8, i.e., eight nontrivial evaluations for exact integration of 
a 15-degree polynomial. 
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The actual evaluation of qj^ thus becomes 

N-1 

^ (Ck - Ck_-|) ^ H- • h(Xj.) sin [q(Xj)] (A-11) 

where 

(A-12) 



The value of N to be used should be consistent with the accuracy re- 
quired and the costs of computation 

To test this procedure, we considered the special case where 





h(x) 


= (1/tt-x)'^ exp[-l/(l/7T-x)] 


(A-13) 




q(x) 


= ttX/(1/tt-x)] 




(A-14) 




a = 


0; b = 1/^ 




(A-15) 


We easily find 












C = 
n 


n/ ^(n+l ) 




(A-16) 


and using N = 


6, we obtain the following 


evaluation 




k 




^k 






0 


0 


- 


- 




1 


1/2tt 


.0225407 


.0225407 




2 


2/37t 


-.0009741 


.021 5666 




3 


3/4it 


.0000421 


.0216087 




4 


4/5it 


-.0000018 


.0216069 




5 


5/6ti 


.0000001 


.0216070 




In this 


particular 


case, the integral can be evaluated exactly. 


We 



have 1 /tt 

I = ir^ / (1-ttx)'^ exp[-Tr/(l-irx)] sin [Tr^x/(l-iTx)]dx 
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Let 1 - Tix = Tiy, and we find 
1 / TT 

r -2 

I = ! y exp(-l/y)sin(l/y-Tr)dy 

"0 

= - exp(-l/y)sin(l/y)dy 

Now let z = 1/y, and we have 

r 

I = - e"^ sin z dz = l/2e'' = .0216070 

^ IT 

and the adequacy of the evaluation is evident. 

We now consider effective means of evaluating the numbers which 

satisfy A-7 and will be bunched more and more closely together as rv>“. 

It may be anticipated that difficulties may be encountered when using 

customary methods of root evaluation. In the applications with which 

this monograph is concerned, g(x) ^ as x -> b” because of the factor 
- 1/2 

(1 - x/b) . If we assume that in the interval (a,b) we have 

[g(x)]2 = ^m(x) 



where m(x) is "approximately constant," we are led to the iterative 
formul a 



bm (cj - a) + a(b - cj ) [g(C^)]^ 
n, (cj - a) + (b - cj) [g(cj)]^ 



where is the jth iterate for C^. One may conveniently take 

c’ = c" , = c , 

n n-1 n-1 

and 

c] = ^ (a + b) 

However, we have not tested this procedure by actual computations. 

For the case of two poles and one zero, explicit evaluations for C 
are easily obtained, as is shown in Appendix B hereof. 
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APPENDIX B 



The Case of Two Zeros and One Pole 



Consider the special case where 

m^(s) = k^(s+a^ )(s+a 2 )/(s+b^ ) 

2 

where k > 0 and 0 < a^ < b-| < a 2 « We have 

y(0) = Pq = k/ a^a2/b^ 

and the inversion is 

F(t) = e'^o - - 

TT 



1 -X t 

^ e sin y dx 



- I |-e"^^siny dx 

TT J X 
^2 



In this case we are able to obtain an explicit representation of the 
values introduced in Appendix A. In both regions of integration we 
have 

y(x) = k/ (x-a^ )(a 2 -x)/(b^-x) 
so that the condition 

w(C^) = niT 

(cf. equation (A-6)) leads to 

^n " ^ [^^1 ^ ^2 ^n^ - ’^(a^ + a2 + - 4(a^a2+b^Yn) 

where 

Yn " n^TT^/k 

The lower sign (-) is to be used for the first integral (from 
a^ to b^ ) and the upper sign (+) is to be used for the second integral 
(from a 2 to «>) . 



(B-1) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 
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APPENDIX C 



Some Additional Inversions 



We have also obtained the following inversions which are not 
related to those considered in the body hereof and which are not listed 
in Reference [3l. They are reported here for record purposes. 

f(s) F(t) 

CO 

cosh/ as e*"\osh/ ac 2 V (-1 )*^cos(nTr/a/b ,-n^Tr^t/b 

— ^ — ' ~ ” ^ / 2 2 ^ 

(s-c) /T sinh/"ET /T" sinh/ be /d c + n n /b 

a/s (n-l)/2 ^ 

S "e (t/a) I^_^(2/TT) (n=l) 

Additionally, we note a rather obvious erratum in 7. 2. 1.2, p. 284, of 
Reference [3] where the factor -2s preceding the sign of summation should 
be -2tt. 
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